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Summary 

An analysis was performed utilizing the model of a wall jet 
for obtaining equations that will predict slot film-cooling 
efficiency under conditions of variable turbulence intensity, 
flow, and temperature. The analysis, in addition to assessing 
the effects of the above variables, makes a distinction between 
an initial region and a fully developed region. Such a 
distinction is important in determining the role that the 
turbulence intensity of the coolant plays in effecting film- 
cooling effectiveness in the area of the slot exit. The results 
of the analysis were used in the correlation of the results of 
a well-designed film-cooling experiment. The result of the 
analysis and experiment was equations that predicted film- 
cooling efficiency within ±4 percent average deviation for 
lateral free-stream turbulence intensities up to 24 percent and 
blowing rates up to 1.9. These equations should be useful in 
determining the optimum quantity of cooling air required for 
protecting the wall of a combustor. 

Introduction 

It is important that the quantity of air used to film cool the 
combustion chamber of a jet engine be sufficient to maintain 
acceptable wall temperatures without affecting engine 
performance more than is necessary. Therefore, an important 
requirement in the design of a high performance gas turbine 
combustor is a good prediction of film-cooling requirements. 
Slot film-cooling prediction approaches applicable to a jet 
combustion chamber must include, in addition to the effects 
of gas flow, temperature, and slot configuration, the effect 
that the turbulence level of the combustion gases and the 
coolant has on the film-cooling efficiency. 

Some work has been done to assess the effect of turbulence 
on the efficiency of slot film cooling. Carlson and Talmor 
(ref. 1), present correlations for conditions of variable free- 
stream turbulence and acceleration. In the analysis that was 
developed to obtain a correlating equation, it is assumed that 
the amount of hot free-stream gas mixed with the film coolant 
is proportional to the mass flow in the boundary layer without 
mass addition. From their correlations of the experimental data 
they obtain values of a mixing coefficient, which are found 
to be proportional to the turbulence level. They determined 
that with higher turbulence intensity of the free-stream the 
higher the mixing coefficient and the larger the quantity of 


hot gas mixed with the film coolant. The correlations of 
reference 1 are for specific values of the turbulence intensity. 
No general correlations in terms of the turbulence intensity 
are presented. 

Juhasz and Marek (ref. 2) performed liner film cooling 
experiments for several slot configurations for conditions in 
the gas turbine combustor of 1 atm and temperatures up to 
1367 K. They found that the film-mass flow rate and the hot- 
stream turbulence level were more important parameters than 
slot geometry. The authors of reference 2 attempted to 
correlate their experimental data by a number of correlating 
equations found in the literature. The correlation that they 
recommend is based on a simple one-dimensional mass mixing 
model where the mixing coefficient is assumed equal to the 
free-stream turbulence intensity. Using this approach, they 
were able to correlate this experimental film-cooling efficiency 
to within ± 20 percent. The results of reference 2 were for 
an average value of the combustion chamber turbulence level 
of 15 percent. The simple mixing model of reference 2 was 
used in reference 3 in an attempt to correlate the film-cooling 
experimental data obtained at four levels of free-stream 
turbulence. In an effort to remove some of the complexities 
of a combustion chamber, the experiment of reference 3 was 
run in a rectangular chamber having a controlled turbulence 
level with a constant free-stream flow and temperature (590 
K). They found that film-cooling effectiveness decreased as 
much as 50 percent for an increase in free-stream axial 
turbulence intensity from 7 to 35 percent. However, they were 
not able to obtain a good correlation between turbulence 
intensity and the calculated turbulent mixing coefficient. The 
authors of reference 3 utilized a constant value of the mixing 
coefficient to predict the experimental data within ±30 
percent. 

The work of Sturgess, et al (refs. 4 and 5), makes significant 
progress in understanding of the role that turbulence plays in 
slot film cooling and in the presentation of correlating 
equations. These correlating equations were developed by 
using functional modifications to an analysis (ref. 6) that 
assumed boundary layer growth of the wall coolant jet. With 
this approach, correlating equations are presented for the data 
of reference 3 (ref. 4). An important feature of these 
correlations is consideration for the potential core region, 
which exists in the region of the slot exit. The assumption is 
made that the film-cooling efficiency is one (rj = 1.0) in the 
potential core region. 

The objective of the present investigation is to extend the 
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initial film-cooling modeling efforts of Juhasz and Marek 
(ref. 2) and Marek and Tacina (ref. 3) into a more 
comprehensive description of the interaction of a coolant jet 
with a high temperature turbulent free stream. By utilizing a 
mixing model, which takes into account some of the features 
of a developing turbulent wall jet, a correlation is developed 
for a prediction of slot film cooling for conditions of variable 
flow, variable temperature, and variable turbulence level. 

A model of a developing turbulent wall jet caused by the 
interaction of coolant and free stream is presented in the section 
Analysis. The model accounts for an initial developing region 
and a fully developed region. A heat and mass balance is 
performed for both regions. The entrance region is divided 
into a fully mixed zone and partially mixed zone to account 
for the effect of the initial turbulence level of the coolant jet. 
A jet model is utilized which grows in proportion to the 
transverse velocity fluctation and the blowing ratio M. An 
analysis of Abramovich (ref. 7) is utilized to determine the 
transition point between the initial and fully developed regions. 
Constants in the developed equations are determined from the 
controlled experimential data of reference 3. Finally, a test 
of the developed correlations is made by comparing predictions 
with the experimental data of reference 3. 

Symbols 

a 0 constant, equation (14) 

b wall-jet thickness 

C M turbulent mixing coefficient, equation (10) 

C 0 constant, equation (7) 

F{R) jet interface position function, equation (25) 
f(x) fraction of coolant mass flow rate, equation (A31) 
g{x) fraction of diffused mass flow rate, equation (A32) 

I v transverse turbulance intensity, based on free-stream 
velocity, equation (18) 

I* transverse turbulence intensity, equation (A50) 

M blowing parameter, p^Jp^U^ 

m mass flow rate per unit span 

N constant, equation (20) 

R velocity ratio 

S injection slot height 

T absolute temperature 

U* jet characteristic velocity 

U gas velocity 

v' fluctuating velocity normal to surface being cooled 

(rms value) 

x distance downstream of slot exit 

wall-jet/coolant interface position, (see fig. 1) 
y 2 wall-jet/free-stream interface position, (see fig. 1) 


(3 initial region factor, equation (4) 
r) film-cooling effectiveness 

V ' film-cooling effectiveness of fully mixed initial region, 
equation (A42) 

6 dimensionless temperature ratio, equation (5) 
p gas density 

0 dimensionless wall-jet thickness, b/S 

\p dimensionless temperature ratio, equation (6) 

co dimensionless flow-temperature grouping, equation 20 

Subscripts 

aw adiabatic wall 

exp experimental 

s coolant 

T total 

0 position of coolant injection, (see fig. 1) 

1 point of transition between initial and fully developed 

region, (see fig. 1) 

I zone I of jet, (see fig. 1) 

II zone II of jet, (see fig. 1) 

°° hot gas 

Superscripts 
— average quantity 

Analysis 

Figure 1 depicts the situation of a secondary fluid (coolant) 
at a low temperature (T s ) being injected at the wall of a 
chamber with a rectangular cross section through which flows 
a high temperature gas. The model of figure 1 simplifies the 
situtation occurring at the liner of a jet engine combustor. The 
effect of the hot and cold streams meeting at the slot exit is 
the creation of a semicontained turbulent jet whose thickness 
b grows with distance x, and whose properties are determined 
by the amount of free-stream fluid entrained m The 
turbulent jet has an initial developing region, sometimes called 
the “potential core’’ region. This region has very efficient 
cooling as indicated by the use of temperature sensitive paints 
on a combustor liner (ref. 8). The length of this region x l is 
assumed for the present analysis to be the distance required 
for the inner edge of the growing/mixing zone I (fig. 1) to 
reach the wall. Since zone II of the initial region consists 
mostly of coolant fluid, some researchers have assumed that 
the wall temperature in this region is at or near the inlet coolant 
temperature T s . Realistically, some fluid (mj) from the 
turbulent mixing zone I will be entrained into zone II and raise 
the temperature in this zone above the value of the inlet coolant 
temperature T s . Since the analysis is for an adiabatic wall, the 
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Adiabatic 

wail 

(a) 




(a) Overall flow model. 

(b) Mass balance control volume. 

Figure 1.— Slot film-cooling flow model. 


wall surface temperature T aw is assumed equal to the mean 
fluid temperature T. In the initial region, the adiabatic wall 
temperature is assumed equal to the mean fluid temperature 
in zone II 

T aw = Tn (x < * 1 ) (1) 


With the above assumptions, a differential mass flow rate 
and energy balance are made of zones I and II (fig. lb), in 
order to relate the film-cooling efficiency to the growth of the 
wall jet. Details of the analysis are given in appendix A. The 
equations derived for the initial region and the fully developed 
region are as follows: 


In the developed region the adiabatic wall temperature is Initial Region x < Xj 

assumed equal to the mixing temperature of the turbulent jet , n (he , nitm| region> the fllm _ cooling effectiveness is 


T aw =Ti (x > Xj) (2) 

The following additional assumptions are made: 

(1) The axial pressure drop of the wall jet is negligible 

(2) The change in the specific heat with temperature is small 

(3) The composition of the mainstream and coolant is the 
same 

(4) The hot gas temperature T '<*, is constant 

(5) The hot stream thermal radiation is negligible 

(6) Interfaces (0-1 and 0-2, fig. 1) of the developing wall 
jet are linear 



1 +0 

/ x\ / (<!>- l)S 
\MS/\ x 


(3) 


Here (3 may be expressed in terms of dimensionless 
temperature ratios 


0exp 


X 

x\ M 


+ 1.0 



(4) 
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0 = 


T — T 

1 aw 1 s 

T - T 

law x =*, l ° W 


0 = 


T — T 

1 aw x = x i iaw 

T — T 

1 OO 1 s 


( 5 ) 


( 6 ) 


or in terms of a dimensionless temperature ratio and turbulent 
diffusion from zone I to zone II 


The efficiency equations (3) and (8) for the initial and 
developed regions may be expressed in a general equation 


where 


1 


V = 


1 + 1— ci 
ms; 


(ii) 


(3 = 


*(+LZl + t.0 


Xi\ M 


x i ~ x \ _ * I ' 

' C 0 K,s ~ — - — In 
x\ ' 



( 7 ) 


Cm — 


Cm — 


x 

MS 


1 + (3 


( 12 ) 


The equation (3 (eq. (7)) is expressed in terms of the turbulent 
mixing coefficient C M (eq. (10)) 


13 = 


„ . x \ x 
CmI — I H — - 


MS 


Xl 


Equation 7 may be used, as shown later, in conjunction with 
equation 3 for predictive purposes. The constant C 0 can be 
determined by using experimental data. 

Developed Region x > Xj 

In the developed region the value of jS is zero and equation 3 
reduces to 


1 


l + ( 

MS 


( 8 ) 


Equation (8) can be expressed in terms of a turbulent mixing 
coefficient C M (refs. 2 and 3) as follows: 


V = 


( 9 ) 


1 Cm 

MS 


where the turbulent mixing coefficient in the present analysis 
is a function of the jet boundary layer growth rate 


X\ — x 


C / 

v -' o l v,s 


In 


X = X\ 


1 - 


X \ , 


with the following conditions 
x > Xi 


(3 = 0 
Cm ~ C M 


Turbulent Mixing Coefficient C 


M 


(13) 


Abramovich (ref. 7, pg. 37) assumes the growth of the jet 
boundary layer b to be proportional to the perturbation 
component of the transverse velocity. This may be expressed as 


db 

dx 


a 0 v 

U* 


(14) 


While both the transverse velocity fluctuation v' and the 
characteristic velocity U* are both likely to vary with x, the 
experimental evidence for turbulent jets suggest a near linear 
jet growth rate. Therefore, it is assumed that edge 0-2 (fig. 1) 
is linear in x. Equation (14) may then be expressed as 


Cm — 


(0 ~ 1 )S 


( 10 ) 


db _ / db\ a 0 v' 0 

dx / U* 

\ / x=0 ° 


(15) 
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Following the approach of reference 9 for the characteristic 
velocity U* of the mixing of two streams with different 
velocities and densities, the characteristic velocity atx = 0 for 
the model of figure 1 is 

U* = + P ^ s - (16) 

Poo Ps 


reference 3, they considered the turbulence of the free stream 
and of the wall jet. Therefore, the effective value of lateral 
turbulent fluctuation v'„ to be used in equation (18) is some 
function of the fluctuation of the free stream >0 and of the 
wall jet v' Q . Based on the recommendation of reference 10 
the following relationship is utilized: 

v 0 = K,o + 0A(\V s [ o ~ V >0 |) (22) 


which can be written as 


or in terms of the defined intensity (eq. (18)) 


* tui + m 
U> 



(17) 


Defining a transverse turbulent intensity as follows: 



(18) 


and substituting equations (17) and (18) into equation (15) 
results in 


db 

dx 


a 0 I v [ 1 + 


(1 +M) 


(19) 


A = ft, « + 0.4( ft, (23) 

with the condition that the value of / v is not greater than both 

■fv.OO tUld /yj. 

Equation (23) suggests that at high free-stream turbulence 
the turbulence of the free stream has the dominant effect on 
the film cooling efficiency than does the turbulence of the wall 
jet. 

Initial Region Length X] 

To determine the length of the initial region x l it is 
necessary to know the position of the interface 0-1 with respect 
to the lip of the coolant slot y x (fig. 1). Abramovich (ref. 7, 
pg. 163), using a momentum and mass balance for the two 
streams, determined analytically the relationship between 
interface positions 0-2 and 0-1. His result can be expressed 
as follows: 


A more general form of equation (19) is 

to aJJ^ 
dx S 


( 20 ) 


where 

oi = (1 + TJT s )/( 1 + M) 

Intergration of equation (20) results in an expression for the 
turbulent mixing coefficient 


^ _ (<£ 1)^ , j AJ y 

0-0 {IyCO) 


( 21 ) 


?2 1 
— 1.0 

— (0.416 + 0.134/?) 

Ps 


(24) 




Pi _ T s 
Ps 

using equation (A25) equation (24) is expressed as 


The work of Juhasz and Marek (ref. 2) assumed the mixing 
coefficient to be proportional to the transverse turbulent 
intensity, but they did not consider the effect of temperature 
and blowing rate indicated in equation (21). 

Ko and Liu (ref. 10) state that the turbulent intensity 
distribution for the main stream and the wall jet have a great 
influence on the film-cooling efficiency. In their calculation 
of the mixing coefficient, according to that proposed in 


Zl 

^2 


= F(R) 


(25) 


- 1.0 


(0.416 + 0.134 R) 


From the geometry of figure 1 the following is obtained: 
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J 2 = 5(0 - 1) (26) 

5 

x l =—x (27) 

y l 

Substituting equations (26) and (27) into equation (25) and 
making use of equation (10) results in the following expression 
for the entrance length: 
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Use of equation (28) requires an initial estimate of jq so that 
a value of the wall temperature at x = jq (T aw ^_ x j may be 
determined. This permits use of equation (28) to calculate jq 
and an iteration routine for the exact value of jq (appendix B). 

Experimental Data 

To enable the equations developed in the section Analysis 
to be suitable for predicting film-cooling efficiencies it is 
required that the constants a 0 , N, and C 0 in equations (21) and 
(13) be determined. This was accomplished by the use of the 
experimental film-cooling data of reference 3. Reference 3 
reported film-cooling data at four levels of free-stream 
turbulence for a constant free-stream velocity (62 m/sec), 
constant free-stream temperature (590 K), constant pressure 
(1 atm), ambient inlet temperature of coolant and four film- 
cooling flow rates. The experiment was designed so that a more 
controlled and simplified version of a jet engine combustor 
would be possible without the complications of the large 
temperature and velocity gradients which exist in an actual 
combustor. In the experiment of reference 3 preheated gases 
were passed through turbulence generators and then through 
a 10.5- by 37-cm rectangular duct, which contained provision 
for slot film cooling on one surface. This surface was insulated 
to permit an adiabatic wall condition. Wall surface 
temperatures were taken of the surface being filmed cooled 
at nine positions downstream of the exit slot position. The 
temperature measurement at the end (9th) location was not 
utilized in this paper because of a wall conduction error. Axial 
and transverse turbulence of the free stream and coolant were 
measured with a hot wire. 

Results 

Constants a 0 and N 

As part of the procedure to determine the constants a 0 and 
N in equation (21) values of the dimensionless jet boundary 



Downstream distance x, cm 


(a) l V OD = 0.026. 

(b) / v c = 0.08. 

(c) =0.13. 

Figure 2.— Calculated dimensionless wall-jet thickness as a function of distance 
from the slot exit for a given free-stream turbulence. 
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(d) / v 00 = 0.244. 
Figure 2.— Concluded. 


thickness (</>) were determined by using equations (3), (8), and 
(28) and the experimental values of reference 3. Calculated 
values of the initial length x { and the dimensionless jet 
boundary layer are given in table I and figure 2. 

Figure 2 indicates that for the lowest three turbulence levels 
there is a linear relationship between the jet boundary layer 
and the distance x. Such linearity is consistent with other 
measurements and analyses of a wall jet. If the initial region 
effect had been neglected and the equation for the fully 
developed region had only been used (eq. (8)) then a nonlinear 
variation in the jet boundary layer would have resulted. In 
terms of the film-cooling efficiency (eq. (3)), the effect of j3 
(eq. (4)) in the initial region is about 10 percent in the 
experiments of reference 3. For the highest value of free- 
stream turbulence (7 V>00 - 0.244, fig. 2(d)) there is a noted 
change in slope at each value of the blowing factor M for values 
of the distance x greater than the distance of the initial region 
jq. Figure 2(d) shows that this slope change occurs at the 
beginning of the fully developed region x = Xj . The slopes 
of the curves of figure 2 for all values of the free-stream 
turbulence varies with the blowing factor. Slopes of the curves 
of figure 2 and values of o> (eq. (20)) are given in table II. 
A slight adjustment was made in the determination of the slope 


for I v oo - 0.244 and M = 0.5 (fig. 2(d)) by having the line 
for the hilly developed region pass through the transition point 
(x = x T ). This introduced a very slight error and made it 
possible to present a single equation for the mixing coefficient 
CM of the fully developed region at the highest turbulent level 
(/v.c. = 0.244). 

Values of the effective lateral turbulence 7 V are also shown 
in table II. These were calculated by utilizing equation (23), 
the experimental values of the free-stream lateral turbulence 
intensity 7 V 00 and the value of the jet lateral turbulence 
intensity I v>s . Average values of the jet turbulence were 
obtained (personal communication with C.J. Marek, Oct., 
1985) of the experiment being utilized in this report at the 
highest blowing rate ( M = 1.85). Values of the jet turbulence 
at other blowing rates were calculated by assuming a constant 
turbulence intensity at the jet exit. 

Equation (20) suggests that a log-log plot of the slopes of 
figure 2 ( d(f)ldx of table II) compared to the product of the 
effective turbulence 7 V and co should result in a linear curve 
of slope N. This is demonstrated in figure 3. The slope of figure 
3 is 0.65 or N = 0.65. With this value for N, a plot is made 
as shown in figure 4 to determine the values of a 0 . 

As a result of the above procedure for determining a 0 and 
N, the following expressions for the mixing coefficients were 
obtained: 

For 

7 V a, = 0.026 
= 0.08 
= 0.13 


TABLE I. -LENGTH OF THE INITIAL 
REGION x x 





X\ /MS 

i 

Calcu- 

lated 


EH 


9.79 

34.44 

39 



23.80 

41.89 

>33 

mm 

1.37 

26.53 

Esa 

>22 

u 

1.85 

35.08 

■ 

>16 

0.08 


5.45 

19.30 

36 



16.35 

29.52 



1.37 

19.97 

23.14 

>22 


1.83 

28.08 

24.34 

> 16 


0.41 

3.50 

13.52 

37 


.86 

11.77 

21.68 

25 


1.36 

16.52 

19.34 



1.86 

20.62 

17.59 

>16 


0.50 

4.00 

B39 

20 


.98 

8.34 


13 


1.40 

11.84 

13.45 

8 


1.84 

14.29 

12.36 

10 


1 














TABLE II. -CALCULATED WALL-JET 
GROWTH RATES 



the equation of the turbulent mixing coefficent for the initial 
and fully developed regions is 

C M = 0.177 (7 v co) 0 - 65 (29) 

and for / v>0o = 0.244 the mixing coefficient for the initial 
region is given by equation (29), while for the fully developed 
region the value of the coefficient is given by 



Dimensionless group of flow, temperature, 
and turbulence, l y u 


Figure 3.— Determination of the constant N (eq. (20)). 



0 .1 o 2 .3 .4 .5 .6 

Dimensionless group of flow, temperature, 
and turbulence, U v u) c - te 

Figure 4.— Determination of the constant a 0 (eq. (20)). 

error in the region of the slot entrance for the data of reference 
3 ranged from less than 1 percent to 4 percent error in the 
wall temperature. The conduction error effect on the 
calculation of C 0 was negligible. To meet the requirement of 
a true constant for C 0 it was necessary to simplify equation 
(13) as follows: 




Constant C 0 

The constant C 0 in equation (13) may be determined by 
equating the value of (3, in equation (13), to the experimental 
value of /3 (eq. (4)), by solving for C 0 , and by using the 
experimental data of reference 3. Estimates of the conduction 
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1.0 


The average value of the constant C 0 determined was 
C 0 — 13.7. The value of the average intensity I*,s °f 
coolant used in calculating C 0 (as determined from the 
personal communication with C.J. Marek, Oct., 1985) was 
0.089. 


Prediction of Film-Cooling Effectiveness 

The following equations are the result of the above analytical 
and experimental considerations: 


F(R) = 


1.0 



(0.416 + 0.134 R) 


1.0 


R = 


U s 


(37) 


1 


V = 


1 + 



(32) 


Fully Developed Region x > X \. — The value of /3 to predict 
the film-cooling efficiency for the fully developed region is 


0 = 0 


(38) 


for 7 vco =0.026 -0.13 


0 



1 +(3 


C M = 0.177 (7» 


0.65 


(39) 


and for 7 V >00 = 0.244 


(33) 


C.M ~ 


0.2451 1 +0.177 X ‘ 


(M 0 - 65 (40) 


Initial Region x < The values of (3 and C M for 
prediction of film-cooling efficiency for the initial region are 
given by 



(1 + M) 



I v = /„,.+ 0.4 (|7 V .„ — A..,|) (42) 

with the condition that 7 V is not greater than both 7 V>00 and I vs . 

With the above equations and the calculation scheme of 
appendix B, predictions were made of the film-cooling 
efficiencies for the flow, temperatures, and turbulence 
conditions reported in reference 3. The resulting predictions 
for film-cooling efficiencies and their comparisons with the 
experimental data of reference 3 are shown in figure 5. 
Comparisons between prediction and experiment are generally 
quite good with an average deviation between prediction and 
experiment being ±4 percent. 


0 - ( T aw x=xi ~ Taw)/ ~ T s ) 


Discussion 


Cm = 0.177 (Iyu) 065 


1 


*i = 



(35) The important role played by the transverse velocity 
fluctuation, as postulated in reference 3, in determining the 
turbulent mixing coefficient C M , is given further confirmation 
in the present work. In addition, the use of a jet model has 

(36) brought out two other variables of importance; namely, the 
blowing ratio M and a temperature ratio TJT S . The forms of 
the turbulent mixing coefficient presented herein do a good 
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Figure 5.— Comparison of predicted film-cooling efficiencies with experiment. 


job of correlating the experimental mixing coefficients of 
reference 3. The developed equations and experimental data 
(fig. 5 ) bring out the important roles played by both the free- 
stream turbulence and the blowing ratio. At low blowing ratio 
M — 0.45 (fig. 5(a)), the effect of turbulence is more 
pronounced than at the higher blowing ratios (e.g. M = 1.9, 
fig. 5(d)). These competing effects need to be considered in 
any design excercises for the optimum cooling of a surface. 

The calculated length of the initial region ^ was found to 
increase with increases in blowing ratio and decrease with 
increases in free-stream turbulence (table I). These lengths are 
significant in their magnitude and for some conditions are 
larger than the length of the experimental measurements 
(ref. 3) utilized in the present work. An estimate of the 
experimental value of may be made by using equations 
(32), ( 33 ), and (40) as a guide for presenting the experimental 


data. For the turbulence range of 7 v oo = 0.026 to 0. 13, a plot 
of I /77 — 1 compared with x/MS should result in the fully 
developed region having a linear correlation that intersects the 
origin. This is shown in figure 6 (a) for the data of reference 3. 
The point at which the fully developed region begins should 
establish the length of the initial region. For the highest 
turbulence level (/ v00 ~ 0.244), equation (40) suggests that 
the linear part of the correlation should intersect th ex/ MS axis. 
This is demonstrated in figure 6 (b). The experimental values 
of x x IMS are given in table I. Except for the lowest blowing 
factor, there is general agreement between the calculated and 
experimental values of the initial region distance parameter 
(X]/MS). 

The correlating equations developed in this paper make a 
distinction between an initial and a developed region. Such 
a distinction, in addition to providing a more accurate jet 
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model, gives some insight into the effect of the initial 
turbulence level of the coolant on film-cooling efficiency. 
Results of the calculations for two coolant turbulence levels 
(I* s ) at high and low free-stream turbulence are shown in 
figure 7. In figure 7, the effect of coolant turbulence on film- 
cooling efficiency is evident in the initial region. The initial 
region for the higher free-stream turbulence (0.13) ends at the 
intersection of the two curves for / voo = 0.13. In the case of 
the lower free-stream turbulence, transition to the fully- 
developed region occurs at a value of x/ MS greater than shown 
in figure 7. Figure 7 is a confirmation of the well known effect 
of slot turbulence on the film-cooling efficiency and suggests 
that the effect is more dominant (overall) for conditions of low 
free-stream turbulence as indicated in references 2 and 6. 



(a) I va> = 0.026, M = 0.45. 

(b) l’ x = 0.244, M = 0.98. 

Figure 6.— Experimental determination of the initial region distance parameter 
X\/MS. 



Figure 7.— Effect of coolant turbulence on film-cooling efficiency, M = 0.9. 

Kacker and Whitelaw (ref. 11) reported no effect of the slot 
turbulence level on the film-cooling efficiency for a 100 
percent variation in slot turbulence. The results of reference 
11 were for an x/S of 50 and 100 which could have been in 
or near the fully developed region where slot turbulence has 
little effect. The experimental data of Sturgess (ref. 5) does 
show an appreciable degree of film-cooling efficiency 
degradation as the coolant turbulence level increases. It can 
be expected that a slot design which produces as low a coolant 
turbulent intensity as possible and which minimizes mixing 
of the two streams should have good film-cooling efficiency 
in the region of the slot exit. In addition, analytical models 
which do not make a distinction between the initial and fully 
developed regions and assume a fully mixed, fully developed 
region for the entire length of the surface to be protected can 
be expected to suffer in validity, especially near the point of 
coolant injection. 

The analytical relations developed in this paper are indepen- 
dent of pressure. This is consistant with the experimental 
findings of Marek (ref. 12) who performed slot film-cooling 
experiments in a high-pressure combustor. Reference 12 
concludes that for the range of pressures investigated (1 to 
20 atm), there is no significant change in the film-cooling 
effectiveness with pressure. 

This paper has emphasized the effect of turbulence intensity 
on jet mixing as it applies to film cooling. The authors of 
reference 3 indicated that entrainment of the hot free stream 
by the coolant jet should also depend on the scale of turbulence 
of the free stream. The work of Oh and Bushwell (ref. 13) 
indicates that in addition to velocity fluctuations, the length 
scale of turbulence also plays an important role in determining 
the turbulent mixing rate of free-shear layers. According to 
the results of reference 13, an increase in turbulent length scale 
(or eddy size) should result in an increased jet mixing with 
a corresponding decrease in film-cooling efficiency. 

Concluding Remarks 

Equations have been developed for predicting slot film- 
cooling efficiency as a function of flow rate, turbulence level, 
and temperature of the coolant and the hot main stream. The 
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analytical results make a distinction between an entrance region 
where there is a potential for efficient film-cooling and the 
fully developed region where complete mixing between the 
hot gases and the coolant is assumed. Constants in the 
developed equations were determined by using the results of 
a well-designed film-cooling experiment with controlled values 
of flow, turbulence, and temperature. The final equations 
based on analysis and experiment were found to predict the 
experimental data within an average deviation of ±4 percent. 


The film-cooling prediction equations presented in this paper 
should be useful in determining the optimum quantity of 
cooling air required for protecting the wall of a combustor 
or of any surface requiring slot film cooling. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 22, 1986 
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Appendix A— Analysis of Film-Cooling Efficiency 


The assumptions used in the analysis are given in the section 
Analysis. Figure 1 depicts the situation to be analyzed, a coolant 
being injected at the wall of a chamber with a rectangular cross 
section through which flows a high temperature gas of the same 
composition as the coolant. The control volume of figure 1(b) 
expresses the differential quantities of mass flow entrained from 
the free stream and the turbulent diffusion of mass from zone 
I (mixing zone) to zone II (potential core). A total mass flow- 
rate balance for zones I and II may be expressed as follows: 


Equation (A7) is solved for the total mass flow rate m Ta of 
zone II by multiplying by the integrating factor 1 l(x x — x) and 
integrating with m Tu = m s at x — 0. The result is 

A + <*'-*> ( X (djA dx (Ag) 

m s \ xj m s } 0 \dxJ (x { - x) 

Assuming a constant dm.il dx and integrating gives 


m T + 


dx 


dm T 

dx = m T - 1 dx 

dx 


which may be integrated to 



m T = m s + Wo, (A2) 

or in terms of zone I and zone II (fig. 1(b)) 

m T = m T] + m T]{ = m s + (A3) 

A differential mass flow-rate balance in zone II (fig. 1(b)) 
results in the following: 


Utilizing figure 1(b), a differential energy balance may be 
applied to zones I and II in the same manner as the above 
differential mass flow balance. The zone II differential energy 
balance is as follows: 

A + A = d_(mXd dx (A1Q) 

dx \Xi — x) dx 


m T\\ 


dmi 

+ — — dx = m Tll + 
dx 11 


dm Tu 

dx 


_ dh 

dx m Tii 

h 11 


(A4) 


integrating equation (A 10) as was done for equation (A7) 
results in 


where dm^dx represents the mass flow rate from zone I to 
zone II due to turbulent diffusion in zone II. The last term 
of equation (A4) accounts for the fraction of the mass flow 
rate in zone II that ends up in zone I as a result of the sloped 
interface 0-1. 

The height h of the interface 0-1 may be expressed as 


fn Tu l n - m a Ts[ 1 ) + (*i “ x ) 

x \j 


' x d(miTi) dx 
J 0 dx ( xi - x) 


(All) 


h = S I 1 - - 
xi 


Therefore, 


(A5) 


The zone I differential energy balance is as follows: 




dx = 


dx 


dx 


(A12) 


dh dx 

h Xi — x 


Combining equations (A10) and (A12) and integrating 
(A6) results in 


with equation (A6), equation (A4) is expressed 


m Tu 7jj + m-j-Ti ~ m s T s -I- (A13) 


dm Tu 

dx 


dx + 


m Tn 

(Xi - x) 


dx = 


dmi 


dx 


dx 


(A7) 


Substituting for in equation (A 12) with equation 

(A 13) gives 
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(m T T\j 

dx 




\xi — x 


m s T s 


X\ — x 


dx = 


dx \x { - x 

d 


dx 


dx 


(A 14) 


yi - — (</>i - i)^ 

*i 


where (j>i = b x /S. 
Since 


Multiplying equation (14) by the integrating factor 
\/(x\ - x) and integrating produces the following: 


m s - U sPs S 

M — U S 6 S / UooPoa 


(A20) 


(A21) 

(A22) 


m T]^I — m sT s 


X 

b mo, — (xj — x) 

X\ 


C x d(miT^ dx 
J 0 dx (x 1 - x) 


use of equations (A19) to (A22) results in the following 
equation for the average temperature in zone I for x < x t : 


(A15) 


m T[ may be determined by utilizing equations (A3) and (A9) 
resulting in 


71 = + 

1+ (^0 1+ M 


+ 


Ul - X ) 


M 


01 ~ 1 


*1 


x ( (01 — 1) 
m, — ( 1 + 


M 


m T = m s b - (^ — x) 

*\ 


( d(mi) dx 
dx (xx — x) 


(A 16) 


f x dmi dx x d(rriiTi ) dx 

J 0 dx (xj -x) J 0 dx (xj -x) 


Substituting for m Tl in equation (A 15) and solving for 7) 
gives 




+ 


i + 


Woo 

X 

m s — 

Xi 


T 

1 CO 


X 

m s — 
Xi 

1 + 

Woo 


(X! - X) 


X 

m s b Wqo 

Xi 


' x dm l dx f x J(wi7j) dx 

, o dx (xj - x) J 0 dx (xj - x) 


(A17) 


(A23) 


Inspection of equation (A23) indicates that for a_ given 
blowing rate M, T s and T x , the variation in x of 7j, will 
depend only on the last term of equation (A23). This last term 
can be shown to be quite small so that 71 is determined by 
the first two terms of equation (A23) which are invariant in 
x. Therefore, a useful simplification for zone I for x < X] is 

71 = constant (A24) 

which permits the following identity in the initial region 

(X < Xj) 


7i 


= T„ 


(A25) 


A mass flow-rate balance of both fluid streams (fig. 1(a)) 
at position x = 0 and at any other position results in 

m T - m s + t/ooP»y 2 (A 18) 

The entrained mass flow rate m* by comparison with 
equation (A3) is 

Wo, = (A 19) 

Because the edges 0-1 and 0-2 (fig. 1(a)) were assumed to 
be linear in x, y 2 has the following relationship: 


Because of the adiabatic wall assumption, the average 
temperature of zone II may be expressed 

Til = T aw (A26) 

Use of equations (A24) to (A26) simplifies equation (All) 
and (A 15) so that together with equations (A9) and (A16), the 
energy and mass flow-rate equations for zones I and II (fig. 1) 
may be expressed for zone I as 

m Tl T aWx=x] = m s f(x)T s + - g{x)m{T aWx ^ (A27) 

m T[ = mj(x) + Woo - g(x)m l (A28) 
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and for zone II as 


(A39) 


mjJ'aw = rn s ( 1 — f(x))T s + gix)m{T aw (A29) 


where 


m Tn = m s { \ -f(x)) + gix) mi 


x 

fix) =- 
X\ 


gix) = 


*1 


1 ) In 


1 -* 


(A30) 


(A31) 


(A32) 




Equation (A30) has the expected form for zone II, since as 
x — Xi, gix) — 0, and m Ta — 0. Combining equations (A27) 
and (A29) and substituting for m Tu and m x (eq. (A3)) results 


in 


m Ty(T aWx=x 


m T\ T s - T a 


1 aw \ ^ T aw — Tq. 

Z - 


ffh 

m T 


(A34) 


Since the film-cooling efficiency rj is defined as 


T - T 

± oo x aw 


- Z 


(A35) 


6 = 


T — T 

1 aw A s 

T — T 

law x=x x law 


Substituting equations (A 18) and (A20) into equation (A38) 
results in 


m Tl _ x 
m s x x 


<h ~ 1 
M 


+ 1 


- i 1 --10 


(A40) 


Since y 2 = b - 5 (fig. 1(a)) equation (A18) may be 
expressed 


m T 


1 


1 + 


MS 


i<t> ~ 1)5 


(A41) 


where (j> = b/S. 

Equation (A36) may be expressed in terms of the efficiency 
r ] ' that would be predicted by neglecting the “potential core’ ’ 
(zone II) (i.e., a fully mixed zone I and II) 




rj =— = — 

m T 


l + m ' T \( Taw x=x l * aw 


m. 


- Z 


(A42) 


where rj — the experimental film-cooling efficiency. 

Using equations (A40) and (A41) results in equation (A42) 
being expressed as follows: 


Equation (A34) becomes 


- Z 


^ m X | f Taw x=x ! 

V m T m T \ - T s 


(A36) 


The last term of equation (A36) is normally neglected in 
film-cooling calculations. This term accounts for the existence 
of “potential core” (zone II). 

Substituting for m Tu in equation (A29), using equation 
(A30), and solving for mi give 


mi — 


gix) 


x 

1 -- 


X\J \Zw x=X{ Zw, 


(A37) 


Use of equations (A30), (A37), and (A3) results in the 
following expression for the total mass flow rate of zone I: 


m Tl = m T - mA 1 ) ( 1 + 0) 


V = 


V 


1 


where 


$exp 


1 + fiexp j_|_^ x \ f i4> ~ 1)5 
MS/ \ X 




(A43) 




M 


X\ 


f (A44) 


d = 


T — T 

* aw i 

T — T 

law x=x x aw 


f = 


T — T 

law x=x 1 Iaw 

T - T 

1 oo X S 


For determining a general value of (3 for equation (A42) 
(A3 8) that will account for turbulent diffusion in zone II, equation 
(A 16) is combined with equations (A 19) and (A20) to obtain 
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0 = 


X ( <t> 1-1 


x,\ M 


+ 1.0 


*1 — X 


m. 


[ x (dm, 

Jo \ dx . 


dx 


(*i -x) 


(A45) 


The turbulent diffusion of mass from zone I to zone II 
0 dm,/dx ) may be expressed in terms of the jet lateral 
fluctuating velocity v' so and the density p, of zone I as 
follows: 


dm i 

T S,oP\ 
dx 


(A46) 


A more exact form of the derivative may be obtained from 
an approximation which utilizes equation (A37) and assumes 
a linear variation in x of the wall temperature T aw 


dm\ _ U s p s S 
dx J X] 

/ pvn A 


(A47) 


exp 


Therefore, the mass flow rate from zone I to II may be 
expressed as 


dm, = C Q v' 0 p { S 
dx x i 


(A48) 


Substituting equation (A48) in equation (A45) results in the 
following: 


0 = 


x /(f), — l 


+ 1.0 


x, \ M 




, 'o*v,s „ 

x, / T„ 


1 - 


X, 


(A49) 


where 


The value of the constant C a may be determined by 
utilizing experimental data. 

Inspection of equations (A44) and (A49) indicates that the 
value of f3 became zero at the end of the initial region 
(x — jcj). The equation for the fully developed region (x > jq) 
is therefore 


1 


V ~ 


1 + 


x\/(ct>- 1)5 


(A51) 


s MS / \ x 

or in terms of a turbulent mixing coefficient C M (refs. 2 and 3) 

1 


V 




(A52) 


where the turbulent mixing coefficient in the present analysis 
is a function of the jet boundary layer growth rate 


Cm — 


(0 - 1)5 


(A53) 


The efficiency equations (A43) and (A51) for the initial and 
developed regions maybe expressed in a general equation 


1 


V = 


1 + [ — ) C M 
MS 


(A54) 


where 


Cm ~ 


Cm — 


x 

MS 


1+0 


(A55) 


with the following conditions: 


JC > JCj 


0 = 0 
Cm = Cm 


j* V S,0 

*v,s 

u . 


(A50) 
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